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Abstract 
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^^ We present a comprehensive analysis of 7Z interference corrections to the weak charge of the 

(~| proton measured in parity-violating electron scattering, including a survey of existing models and 
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a critical analysis of their uncertainties. Constraints from parton distributions in the deep-inelastic 
region, together with new data on parity-violating electron scattering in the resonance region, 
result in significantly smaller uncertainties on the corrections compared to previous estimates. At 



K^ the kinematics of the Qwcak experiment, we determine the 7Z box correction to be KeD!^ = 

^^ (5.60 lb 0.36) X 10^^. The new constraints also allow precise predictions to be made for parity- 

^ . violating deep-inelastic asymmetries on the deuteron. 
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I. INTRODUCTION 

Modern low-energy experiments at the precision frontier provide important alternatives 
to high-energy tests of the Standard Model currently being performed at the Large Hadron 
Collider (for recent reviews, see Refs. P-[3]). One such experiment is the parity- violating 
(PV) elastic electron-proton scattering measurement that was recently carried out by the 
Qweak collaboration at Jefferson Lab [4J, which aims to determine the proton's weak charge 
Q^ to within 4%. At tree level, the weak charge is related to the weak mixing angle, sin^ 9]y, 
by QvK ~ -'- ~ 4sin^6'vy. By scattering low-energy polarized electrons from an unpolarized 
hydrogen target, Qweak measured the asymmetry between the cross sections for right- and 
left-handed electrons, 

Apv = -^— -, 1) 

where ax is the cross section for a right-hand (helicity A = +1) or left-hand (helicity A = — 1) 
electron. At small four-momentum transfer squared t, the asymmetry is related to Q^^ by 

where Gp is the Fermi constant and a is the fine structure constant. Including radiative 
corrections, the proton's weak charge can be written as [6] 

where sin^ Ow{^) is the weak mixing angle at zero momentum, and the correction terms Ap, 
Ag and A^ are well understood and have been computed to sufficient levels of precision [H] . 
Similarly, the work of Refs. ^\^ has established that the electroweak box diagrams Dw/iy 
and \I\zz are known within Qweak uncertainty limits. 

Until recently it was also believed that the interference 'jZ contribution, illustrated in 
Fig. [T| was known to sufficient accuracy for the Qweak experiment. This correction is defined 
in terms of the electroweak amplitudes as |10j 



O,zi0) = Ql. -^ — (, (4) 
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FIG. 1: Interference 7Z box (left) and crossed box (right) diagrams. The wavy and dashed hnes 
represent the exchanged 7 and Z bosons, with the electron, hadron and virtual photon momenta 
labeled by /c, p, and g, respectively. 

where M.^ is the electromagnetic Born amplitude, JVv^, is the parity- violating part of the 
Born Z exchange amplitude, and M\z ^^ ^'^^ parity-violating part of the 7Z interference 
amplitude (including the contributions with the 7 and Z interchanged). A groundbreaking 
contribution was made by Gorchtein and Horowitz [llj, who showed, using a dispersion 
relations approach, that the D^^ term was strongly energy dependent and was much larger 
at Qweak energies (~ 1 GeV) than previous estimates had assumed |[6j. More importantly, 
the uncertainty on this correction was such that it could significantly affect the precision 
aims of the Q„oak measurement. 

Subsequent analyses by Sibirtsev et al. [12J and Rislow and Carlson [13] generally agreed 
with the overall scale of the correction found in Ref. [H], but disputed the magnitude of the 
uncertainties. In a follow-up study, Gorchtein et al. [H] performed a more detailed analysis 
of the model dependence of the D^^ contribution, correcting several errors from the original 
analysis [H] , but still quoted uncertainties twice as large as those in Refs. [121 [13] . 

Since the interpretation of the Qweak results depends on having a sound understanding of 
the d-yz correction, the lack of consensus about the magnitude of its uncertainty is obviously 
problematic. To move beyond this impasse, in this paper we revisit this problem with the 
aim of resolving the disagreements. 

We begin our discussion by outlining in Sec. [IT] the dispersion relation formalism used to 
compute the 7Z corrections in terms of jZ interference structure functions. The latter are 



the main input into the calculations and are reviewed in detail in Sec. Ill In particular, we 
discuss the uncertainties in determining the 7Z structure functions from electromagnetic 
data for both the resonance and nonresonant background contributions. Constraints from 
parton distribution functions in the deep-inelastic scattering (DIS) region and new data from 



the parity-violating electron-deuteron scattering experiment E08-011 at Jefferson Lab [15 



in the resonance region are used in Sec. IV to hmit the uncertainty range in models for the 
7Z structure functions, and to provide more reliable bounds on the box corrections. The 
resulting D^^ correction is presented in Sec. \V\ where we contrast the revised uncertainties 
with those estimated in previous unconstrained analyses. Predictions are also made for 
parity-violating deuteron asymmetries in the deep-inelastic region, as well as for the recently 
completed inelastic measurement by the Qweak collaboration [16]. Finally, we draw some 



general conclusions from this analysis in Sec. VI and explore possibilities to further reduce 



the uncertainties on the jZ corrections in the future. 

II. DISPERSIVE ANALYSIS OF PARITY- VIOLATING ELECTRON-HADRON 
SCATTERING 

The jZ interference correction O^z can be decomposed into two parts, arising from the 
vector electron — axial vector hadron coupling to the Z boson (Df^)) ^^^ from the axial 
vector electron — vector hadron coupling to the Z (D^^^), 

n,ziE) = n^.iE) + n];,iE). (5) 

At very low energies, such as those relevant for atomic parity violation experiments [T71 
\TE\ , the n;^2term dominates, while the contribution from the D^^is negligible. At the 
energy of the Qweak experiment, however, both terms provide significant contributions. The 
n^z corrections were first computed some time ago by Marciano and Sirlin [3 |8] and were 
updated recently within a dispersion relation framework by Blunden et al. ^9\ [20], with 
reduced errors. The vector hadron correction, DY^ , which is subject to significantly larger 
uncertainty, will be the focus of the rest of this analysis. We will consider only the inelastic 
contribution to D^z ] the elastic contribution has previously been considered in Refs. [3 El 

122] and is strongly suppressed by an additional factor Q^. 

For forward scattering, the dispersion relation for the real part of OYz is given by 

2E„ f^ „, 1 



Ken^^(E) = —Vj^ dE'^j^-^^Qmn':;,{E'), (6) 

where V denotes the principal value integral, and we have used the fact that D^^ is odd 
under the interchange E' ^ —E'. From the optical theorem, the imaginary part of the PV 



^Z exchange amplitude can be written as [T0HT2 



where Q^ = —q^ represents the virtuahty of the exchanged boson, and the integration 
variable k' = k — q. The jZ lepton tensor is given by 

L;^f = u{k, A) ((7^7^ - 9ll,l5) r 7. u{k, A), (8) 

where the vector and axial-vector couplings of the electron to the weak current are Qy = 
— (1 — 4sin^^i4/)/2 and g"^ = —1/2, respectively, and A is the lepton helicity. The hadronic 
tensor for a nucleon initial state is defined as 

X 

where Jt^ and J^ are the electromagnetic and weak neutral currents, respectively, and px 
is the four-momentum of the hadronic intermediate state X. Using isospin symmetry, the 
matrix elements of the vector component of the Z current for a proton target can be related 
to the proton and neutron matrix elements of the electromagnetic current by 

(X|J^b) = {1- Asm' ew){X\Ji;\p) - {X\Jt;\n), (10) 

neglecting the small contribution from strange quarks. In general, the hadronic tensor can 
be decomposed in terms of the 'jZ interference structure functions F^ as 

MW^; = -g^^-F^^ + ^^F^^ -le'^'^P^^^Ff, (11) 

^ p ■ q 2p ■ q 

where p is the four-momentum of the target hadron. Note that the structure functions F^ 
and F2 contribute to the vector hadron contribution, while the F^ structure function ap- 



pears only in the axial vector hadron correction. Combining Eqs. (pj) and ( 11 ), the imaginary 
part of the nY^ correction becomes 



^-'^)-(73W/„>^''^/'"^''«^r^f7A? 



X 



z 
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(12) 



p"/Z . S (Qmax ~ Q ) p-yZ 

where s = M^ + 2ME is the total center of mass energy squared, W^ = (M + m^r)^ is the 
mass at the pion threshold, and Qmax = 2ME{1 — W'^ /s). Following Ref. [19], we include 



in Eq. (12) the Q dependence in a[Q ) arising from vacuum polarization contributions. 



The most important inputs into Eq. (12) are the ^Z interference structure functions 



Fl , which are functions of two variables, usually taken to be Q^ and the Bjorken scaling 
variable x = Q'^/2p ■ g, or alternatively Q^ and W"^. Unfortunately, these functions are not 
well determined experimentally. Although there are some data on F^' and F^ at high W 
and Q^, in the low-W^ and Q^ region, which is crucial to the dispersion integrals, there is 
little or no information. Unlike the electromagnetic structure functions, which can be fit to 
the ample data available, the F^ must be expressed through models. Given that it can 
be difficult to resolve the accuracy of the models, the controversy in the literature over the 
3fJe D^^ contribution is not surprising. 

For later reference, we note here that the Fi and F2 structure functions, for either 7Z or 
electromagnetic (77) scattering, can be related to the transverse (ar) and longitudinal (o"l) 
electroweak boson production cross sections as 



Fi{W\Q') = (^^^^^)aT{W\Q% (13a) 

where u = E — E' is the energy transfer. For convenience one often defines the longitudinal 
structure function as the combination of Fi and F2 structure functions given by 

Fl = (1 + ^] F2 - 2xF,, (14) 

where the prefactor can also be written as (1 + 4x^M^/Q^). 

III. 7Z INTERFERENCE STRUCTURE FUNCTIONS 

Most of the uncertainty in the calculation of the D^^ correction arises from the incomplete 
knowledge of the 7Z structure functions. There have been extractions of F^ and xF^ 
from neutral current DIS by the HI collaboration at DESY [23] at very high Q^ (60 < 
Q^ < 50,000 GeV^) and small x (0.0008 < x < 0.65) using longitudinally polarized lepton 
beams at HERA. However, these data have little overlap with the region of most relevance 
for the dispersion integral, which receives contributions primarily from high x and low Q^, 
where there are no direct measurements. Consequently, one must appeal to models of the 
interference structure functions to estimate D^^ . 



In this section we review the models used in the hterature for the 7Z structure functions, 
before presenting our constrained model, which we refer to as the Adelaide- Jefferson Lab- 
Manitoba (AJM) model. The construction of the models involves first choosing appropriate 
electromagnetic structure functions F^^ , and then transforming these to the 7Z case. In 
describing the structure functions, or equivalently the virtual boson-proton cross sections 



(Tt,l in Eqs. (13), it is convenient to separate the full range of kinematics into a resonance 



part and a smooth nonresonant background, 

(res) . (bgd) /^ f.\ 

The CTrp'^l term includes a sum over the prominent low-lying resonances, while Cj^f^ is 



determined phenomenologically by fitting the inclusive scattering data [211 [25]. Although 
such a separation is inherently model dependent, as only the total cross section is physical, 
it nevertheless provides a useful way to parametrize the somewhat different behaviors of the 
cross sections in the low- and high-iy regions. 

For completeness, the following list summarizes the models for the 'jZ structure functions 
that have been discussed in the literature: 



(i) color-dipole model [261 127], referred to as "Model I" in Gorchtein et al. (GHRM) 

(ii) vector meson dominance ( VMD) + Regge model [281 129] , referred to as "Model II" by 
GHRM [H]; 

(iii) Sibirtsev et al. (SBMT) model [12], based on the Regge parametrization of Capella et 
al. [30]; 

(iv) Carlson and Rislow (CR) model [T3l l3T]. 

The models [1214X^1 [31] differ primarily in the treatment of the background contributions 
ctt^I for the jZ interference, the uncertainty on which is the main source of disagreement 
between the various estimates of D^^ . For the resonance region, all of the models (with 
the exception of SBMT [T2]) use the Christy and Bosted (CB) parametrization [21] of the 
electromagnetic structure functions at low W, but differ in how these are transformed to 
the 7Z case. Note, however, in both Model I and Model II of GHRM some of the resonance 
parameters in the CB fit are modified to better match the choice of background contribution 
[14] . In the following we discuss both the resonance and background content of these models 
in more detail. 



A. Resonances 

The CB parametrization [21] of F^l^ fits the resonance region electron-proton scattering 
data in terms of the seven most important resonances (P33(1232), Pii(1440), Di3(1520), 
511(1535), 5*15(1650), Fi5(1680) and an / = 3 state with mass 1934 MeV), and generally 
agrees with the data to within 5%. The CB fit is used as the basis for the resonance models 
of Carlson and Rislow [31j, and Gorchtein et al. |T3], with the latter using slightly modified 
parameters for cr^'^ in their Models I and II. Sibirtsev et al. [12], on the other hand, perform 
their own fit of the data, incorporating the four resonances P33(1232), Di3(1520), Fi5(1680) 
and ^37(1950), and also obtain a reasonably good description of the data. 

Modifying the electromagnetic structure functions to obtain their interference analogs 
involves modifying the contribution from each resonance P by a ratio that takes into ac- 
count the differences between the electromagnetic and weak neutral transition amplitudes. 



according to Eq. (|10|). For the transverse cross section GHRM define this ratio for a proton 
as 



e^ = ^ = (l-4sin2M-2/R, (16) 



^T,R 



where 






AP An* _|_ AP An* 

^R 1 ^R 1 ^ ^R 3 ^R 3 
= ^'-2 ^'1 ^'-2 ^'-2 (17) 



-K, 2 I "-'2 



with A^ -^ the transition amplitude from a proton or neutron to a resonance R with helicity 
A = ^ or |. The amplitudes A^^ are assumed by GHRM to be Q^ independent, and their 
values determined from electromagnetic decays at Q^ = [32]. The ratio for the longitudinal 
cross section is taken to be equal to the transverse ratio in both Models I and II of GHRM. 



Carlson and Rislow [3l] use a similar ratio to that in Eq. (16) (which they label as Cr), 
but include in addition a Q^ dependence in the amplitudes derived from the MAID unitary 
isobar model [33]. For comparison, CR also calculate the transition amplitudes using a 
constituent quark model [13j . 

Finally, Sibirtsev et al. [12] use the conservation of the vector current and isospin symme- 
try to set the ratio for isospin-3/2 states to (1 + Q^) ^ (2 — 4sin^ 9w)- For the isospin-1/2 
resonances, such as the Di3(1520), SU(6) quark model wave functions are used to estimate 
the ratio of couplings. The similarity of the magnitudes of the weak and electromagnetic 
couplings was used by SBMT to justify approximating the ratio C,Rhj 1. 
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B. Background 

1. Electromagnetic structure functions 

Although the CB parametrization ^^ includes a background cr^| at low W (W < 
3 GeV), to describe the nonresonant contributions to the electromagnetic structure functions 
at W > 3 GeV requires a model for the background which is also valid at large W. In the 
calculation of GHRM [H], the color dipole model from Cvetic et al. [26l |27] is used for 
Model I, while the VMD+Regge model of Alwall and Ingelman ^] is employed for Model II. 
Since the latter was shown by GHRM to introduce the largest uncertainty in D^^ , it will 
be the main focus of our attention. 

According to the VMD hypothesis, the interaction of a photon 7 with a hadron proceeds 
through transitions to vector mesons V (with \^ = p, cj or 0), with strength A/47ra//y, where 
fv is the electromagnetic decay constant of V . The three vector mesons saturate around 
80% of the total photoproduction cross section [28j. The remainder is usually attributed to 
contributions from higher masses, which are modeled by a continuum of states starting at 
mass mo ~ 1.4 GeV [28]. (In the case of the color dipole model [261 E3 El], the photon is 
assumed to interact with the hadron through coupling to uncorrelated qq states instead of 
mesons.) Following Ref. [29], we neglect the off-diagonal terms in the mass integral, which 
is known to be a good approximation for scattering from nucleons [33] . The transverse and 
longitudinal virtual photon-nucleon cross sections can then be expressed as 



,VMD 

Tj, — a^N 



^VMD 



\^ 1 1 



V 



;i + gVO' "i + QV"^o 



;i8a) 



Q^m 



+ Qymiy 



. V ^ 



:i8b) 



where a^N is the real photon-nucleon cross section, and the constants ry ~ l/Zy represent 
the relative contributions from the individual vector mesons V, with r^ = 1 — ^ ry being 

V 

the continuum fraction [29]. Phenomenologically, the ry values are determined as ry = 
{0.67, 0.062, 0.059} for V = p, u and 0, respectively [36]. As we shall see below, re plays a 
critical role in determining the uncertainty on the interference cross sections. The parameters 
C,v and C,c allow for different behavior of the transverse and longitudinal components of the 



vector mesons, although in practice these are usually set equal, ^y = ^c, in order to fit the 
available data. Note that despite the apparent l/Q"^ dependence in the second term of aj^^ 



in Eq. (18b), one can verify by expanding the logarithm for small Q^ that the longitudinal 
cross section does in fact vanish in the Q^ — )■ limit. According to Regge theory, the real 



photon cross section can be parametrized as a sum of two terms [37j 



a^N 



^•y S^ ~r -Dj S^ , 



7 -^7 



(19) 



where s^ = W"^, with the exponents e and r] giving the energy dependence of the Pomeron 
and Reggeon terms, which have coefficients A^ and B^, respectively. 

In the model of SBMT, the background is parametrized according to the structure func- 
tion fit of Capella et al. [3D], with several parameters adjusted to better describe recent 
data, as discussed in Ref . [12] . The parametrization of the F^'^ structure function, which is 
valid for all Q^, is again given by a sum of Pomeron (P) and Reggeon R exchange terms. 



f:1\x,q'') = Apx-^(i-xr+^ 



Q' 



nl+A 



Q2 + A2J 



+ A^x'-'^^il 



X] 



Q' 



Q2 + Ai 



"R 



(20) 



where A and n are both functions of Q^, and Ap, Ap, Ar, Ar and or are fit parameters 
The F^'^ structure function is obtained by SBMT from a parametrization of the ratio 



of longitudinal to transverse cross sections. From Eqs. (13) this can be written as 






1 + 



(21) 



Q^ J 2xFi 

which is parametrized by a sum of exponentials |12j . 

While the above models use the same background parametrization over the entire range 
of kinematics, CR [T3| [31] on the other hand, divide their dispersion integral into three 
distinct regions, each described by a different model. In particular, the resonance region at 
low W is described in terms of the CB fit to cr^*^ and aj,f^ [21], while for the high-H^, 
low-Q^ region, CR use the Capella et al. structure function parametrization. For high W 
and high Q^, a partonic description is employed using the CTIO global fit [3H] of parton 
distribution functions (PDFs). 
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2. ^Z structure functions 

To construct the nonresonant background contributions to the transverse and longitudi- 
nal 7Z cross sections, the electromagnetic cross sections need to be rescaled by the ratio 
^tlI'^'tl^ as for the resonance components. For Model II of GHRM [H], a generalization 
of the VMD model is used, assuming the 7Z cross section for vector meson V is given by 
the analogous 77 cross section scaled by the ratio ny of weak and electric charges. 



a. 



T,L 






(22) 



where 



Kp = 2 — 4sin^6'vy, 



Ku 



-4 sin 9w, 



K^ = 3 — 4 sin^ 9 



w 



(23a) 
(23b) 
(23c) 



correspond to the isovector, isoscalar and strange quark components of the electroweak 
current, respectively. This allows the ratio of 'jZ to 77 cross sections to be written as 



a. 



T,L 



77 



a. 



T,L 



I + Rl^\Q^) + Rl^\Q^) + R'd\Q^) 



yT,L ■ 



where Ry is the ratio of cross sections for V and the p meson. 



R 



T,L 
V 



a. 



(T- 



77{V^) 
77(p) 



fl n+ Q^lm; 
f^\l + Q^/ml 



-,T,L 



The corresponding ratio Rq of the continuum to p contributions is given by 



Re 



Re 



Tc fi + QV^]^' 



r. M + QVmg 









° ln(l + Q^/ml) 



QVml 



Q^/ml 



QVmlf 



(24) 



(25) 



(26a) 
(26b) 



with the continuum mass parameter set to itlq = 1.5 GeV [11]. The parameters Kq in 



Eq. (24) denote the ratios of the jZ and 77 continuum contributions to the cross section. 



Unlike for the discrete vector meson terms, the VMD model does not prescribe a simple 
charge ratio factor to modify the continuum part of the cross section. In view of this. 
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GHRM proceed by assigning a 100% uncertainty on this contribution. As we will see below, 
this assumption gives the largest contribution to the uncertainty on D^^ . 

For Model I of GHRM, the same general form for the 7Z cross sections is used as in 



Eq. (24), but with different individual contributions Ry . Whereas in Model II, the Ry 



are functions of Q^, in Model I these become constants with relative strengths determined 
by squares of quark electric charges, with the continuum contribution associated with the 
J/ip meson [14J, {p : u : (p : J ftp} = {1:1/9:2/9: 8/9}. Similarly, a 100% uncertainty is 
assumed for the J/ip term in Model I. 

In the SBMT model [12], the jZ structure functions at low x are approximated by their 
electromagnetic counterparts. This is motivated by the approximate flavor independence of 
sea quark distributions in the low-x region, and the similarity of the sum of the electroweak 
couplings for 3 quark flavors, ( ^ ^g9v) / [ Xlq ^g ) — 2 — 4 sin^ 6w ~ 1 p4], where Cg and gy 
are the electric and weak vector charges of quark q, respectively. At larger x (x > 0.4), how- 
ever, SBMT compute F^ using a ratio of leading twist (LT) structure functions computed 
from the MRST parton distribution functions [55] , 

F^' = ( ^ 1 Fr. (27) 

At these x values, SBMT note that the flavor dependence of the parton distributions renders 
the interference function ~ 30% — 40% smaller than the electromagnetic, which therefore 
provide an upper limit on F^ . 

Finally, for the CR model [131 EI] the method for modifying the 77 background cross 
sections depends on the kinematic region of W and Q^. In the resonance region, CR take the 
average of the high energy (x — )■ 0) limit (u = d = s), in which F^ / F^'^ = 2 — 4sin^ 9w, and 
the SU(6) quark limit {u = 2d, s = 0), in which F^ / F^^ = 5/3 — 4sin^ 9w, to convert the 
electromagnetic background from the CB structure function parametrization [21]. For the 
low-Q^, high-VT region, CR apply the same ratio to the Capella et al. [30] parametrization 
as SBMT, while in the DIS region they compute the F^ structure functions directly from 
LT parton distributions [38] . 

Using these models for the resonance and nonresonant background contributions to the 
7Z structure functions, the analyses of GHRM [2], SBMT jJ2j and CR [13] estimate the 
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n^ - 


= (5.4 ±2.0) X 10^3 


GHRM 


□;;. = 


= (4.7+^i)xl0-^ 


SBMT 


□?■. = 


= (5.7 ±0.9) X 10"^ 


CR] 



D^^ correction at the Qweak energy to be 

(28a) 
(28b) 
(28c) 

respectively. The GHRM result for the central value of 3fJe Dl^^ is the average of Mod- 
els I and II, but with the dominant background error taken from the larger of the two, 
in this case Model II. The GHRM analysis also estimates the effect of the t dependence 
of the n^^ correction, from t = in the dispersion formalism to t = —0.03 GeV^ in the 
Qweak experiment, finding a decrease of approximately 1.3%, with a similar uncertainty on 
the correction at the Qweak point. 

The central values of all the calculations agree within the quoted uncertainties; however, 
the error on the GHRM value is twice as large as those on the SBMT and CR calculations, 
even though the SBMT estimate includes a fairly conservative uncertainty on the input 77 
structure functions. Given the importance of the D^^ correction to the extraction of the 
weak mixing angle from the Qweak measurement, it is vital that the origin of this difference 
be understood, and ways of further reducing the uncertainty explored. 

C. Adelaide Jefferson Lab Manitoba model 

To proceed with our analysis of the 7Z correction, we define here the ingredients of our 
Adelaide-Jefferson Lab-Manitoba (or AJM) model, within which we will study in detail the 
various contributions to 3fJe Dl^^ and their uncertainties. We draw on the valuable experience 
obtained with the existing models pT]4T^l3T] . and incorporate into the AJM model some of 
the more robust features of the previous analyses. Most importantly, we consider additional 
constraints from existing data on some of the model parameters which were unconstrained in 
the earlier work. We will find that indeed data on PDFs near the resonance-DIS transition, 
together with new results on inclusive parity- violating electron scattering asymmetries, place 
significant constraints on the models, in particular on the background contribution. 
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W^ (GeV^) 



FIG. 2: (color online) Kinematic regions contributing to the D^^ integrals in the AJM model: 
Region I (blue shaded) at low W and low Q^ is described by the CB F^l ^t [24j, transformed to 
the jZ case; Region II (red shaded) represents the high-VF, low-Q^ domain as in Ref. |29j (or the 
GHRM Model II [l4j), transformed to "yZ; and Region III (green shaded) at high W and high Q^ 
is described by global PDF fits to high-energy scattering data [40j. 

1. 'j'j structure functions 



Following CR [I3l [3lJ, we divide the integrals in Eq. (12) into distinct regions of W"^ 
and Q^, using specific models to parametrize the ^Z structure functions in each region. 
This is illustrated in Fig. [2| where the W^ and Q^ divisions and the models describing 
them are indicated. Although the boundaries between the regions are clearly defined, the 
models themselves overlap, allowing important checks to be made on the continuity of the 
descriptions across the boundaries. 

For the input 77 structure functions, we use the CB parametrization [24J to describe the 
\ow-W region (Region T) at W-,, < W < 2 GeV for all Q"^ up to 10 GeV^. In fact, the strong 
suppression of the resonance transition form factors with increasing Q^ results in negligible 
resonance contributions already beyond Q^ ^ 2 GeV^. Since the CB fit also describes data 
up to W"^ = 9 GeV^, we use it in the higher-PF region for Q^ < 2.5 GeV^, as indicated by 
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FIG. 3: (color online) Proton F2^ structure function versus W'^ at fixed Q^ = 0.05, 0.5, 1.5 and 
2 GeV^ for the CB fit [21] at low W (blue solid) and VMD+Regge parametrization [29j at high 
W (red dashed). The boundary between these (corresponding to Regions I and II in Fig. ^ is 
indicated by the vertical dashed line at W'^ = 9 GeV^. 

the blue shaded area in Fig. |2} 

At higher W, corresponding to kinematics where Regge theory is apphcable, the 
VMD+Regge model of Alwall and Ingelman [29] is combined with a modified CB resonance 
contribution (c/. Table 11 of Ref. [H]) to describe the structure functions for W^ > 9 GeV^ 
and Q^ < 2.5 GeV^ (Region 11, red shaded area in Fig. ^. Of course, at these values of 



W the resonances will contribute very little to the dispersion integral in Eq. (12), which 



will be contaminated by the background contribution. This model also forms the basis for 
Model 11 of GHRM [H]. The matching of the CB and VMD+Regge parametrizations at 
the boundary between the low-PF and high-VT regions is illustrated in Fig. [3] for the F2"' 
structure function as a function of W"^, at several fixed values of Q^, from Q^ = 0.05 to 
2 GeV^. The agreement between the two models in the region of overlap is clearly excellent. 
For the structure function in the VMD+Regge model, we have assumed a conservative 5% 
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FIG. 4: (color online) Proton ¥2"^ structure function versus W'^ at fixed Q^ = 2.5, 5, 8 and 10 GeV^ 
for the CB fit [21] at low W (blue solid) and the ABMll PDF parametrization [30] at high W 
(green dotted), with the boundary between Regions I and III at W"^ = 4 GeV^ indicated by the 
vertical line. For the Q^ = 2.5 GeV^ panel, the matching with the VMD+Regge model |29j (red 
dashed), corresponding to the boundary between Regions I and II, is indicated by the vertical line 
at W/2 = 9 GeV^. 

uncertainty, similar to that for the CB parametrization. 

In the DIS region at high W and high Q^ (green shaded area in Fig. pi), the structure 
functions can be computed in terms of global PDFs, for which we use the next-to-next-to- 
leading order (NNLO) fit by Alekhin et al. (ABMll) |10]. This fit includes both leading twist 
and higher twist contributions, allowing for descriptions of data for Q^ > 2.5 GeV^ and W > 
1.8 GeV, which overlaps partially with the CB [21] and VMD+Regge [22] parametrizations. 
(Other similar global fits, such as those in Refs. [HHIH], give very similar results, and 
differences between the parametrization generally lie within the PDF uncertainties.) The 
transition between DIS kinematics (Region III) and the models describing the lower-iy 
and Q^ regions is illustrated in Fig. E^ for Fi^'^ at Q^ = 2.5 GeV^ (where the transitions 
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FIG. 5: (color online) Proton Fg^''^ structure function versus Q^ at fixed W'^ = 4, 6, 9 and 12 GeV^ 
for the CB fit [24j (blue solid), the ABM 11 PDF parametrization [40j (green dotted), and the 
VMD+Regge model [29] (red dashed), with the boundaries between Regions I, II and III indicated 
by the vertical lines at fixed Q^. Note that the small disagreement between the VMD+Regge 
model and the PDF parametrization for Q^ = 2.5 GeV^ appears only at larger W'^ values where 
the contribution to the dispersion integral is small. 

between all three parainetrizations is shown at W"^ = 9 GeV^) and at higher Q^ values, up 
to Q^ = 10 GeV^, for the transition between Regions I and III. Again, the models generally 
match very well across these kinematic boundaries. 

The boundaries between the three regions can also be displayed for fixed W^ as a function 
of Q^, as illustrated in Fig. [sj The matching of Regions I and II for W"^ = 4 GeV^ shows 
excellent agreement between the CB |21] and ABMll PDF [40] parametrizations at Q^ = 
2.5 GeV^. At the highest W value at which the CB fit is valid, W'^ = 9 GeV^, the agreement 
between the models describing all three regions is also quite good. For larger W [W"^ > 
10 GeV^) the VMD+Regge model p9j slightly exceeds the PDF parametrization. However, 
this generally occurs at the edge of the kinematic boundary between Regions II and III, 
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where the contribution to the imaginary part of D^^in Eq. (12) is very snialL 



2. ^Z structure functions 

Having detailed the forms of the electromagnetic structure functions, we now turn to 
their 7Z interference analogs. For the low-VT/low-Q^ region dominated by the nucleon 
resonances, the transverse and longitudinal 77 cross sections parametrized in the CB fit 



are modified using the ratio ^r in Eq. (16), with the parameter y^ determined from the 



proton and neutron helicity amplitudes as in Eq. (17). This follows closely the approach 
adopted by GHRM [14J, but, importantly, differs in the way the uncertainties on the helicity 
amplitudes A^y^ are determined. 

In particular, GHRM combined the uncertainties on the amplitudes by adding extremal 
values of each, which implicitly assumes a uniform error distribution rather than the standard 
Gaussian one. Adding errors linearly clearly overestimates the uncertainties, and in the 
AJM analysis we adopt the more conventional Gaussian distribution to add the errors in 
quadrature. (When combining all of the uncertainties on the final 3fJe D^^ value, however, 
GHRM add the errors in quadrature.) In Table ll] the y^ values for the proton and their 
uncertainties computed using both methods are shown for comparison. For completeness, we 
also list the y^ values for the neutron and deuteron, with uncertainties added in quadrature, 
which will be needed in subsequent sections. For the isospin-| P33(1232) and ^37(1950) 
resonances, the uncertainties on the helicity amplitudes are given by the Particle Data 
Group (PDG) |l6] as zero. To be conservative, however, we follow GHRM |13] and include 
a 10% uncertainty on the P33(1232) and a 100% uncertainty on the ^37(1950) resonance 



Note that in Table U and in our numerical calculations we make use of the latest values 
of the helicity amplitudes from the PDG [46J. However, when comparing directly with the 
GHRM analysis [H] we will refer to the earlier, 2010 PDG values [32] that were utilized by 
GHRM for the Di3(1520) and Pii(1440) resonances. The yR ratios using these earlier values 
are listed in parentheses in Table [T| but with errors evaluated using Gaussian distributions. 

For the nonresonant background, the models describing the electromagnetic structure 
functions are transformed to the 7Z case according to the kinematic region considered. For 
the region of low Q^ but high W , the cross section in the VMD+Regge model p9] is modified 
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TABLE I: Electromagnetic to 7Z resonance cross section transformation ratios yn from Eq. (17) 



for the proton, neutron and deuteron in the AJM model, compared with the proton ratio in the 
GHRM model |14j . The AJM model values in parentheses use helicity amplitudes from the earlier 
2010 PDG [32], as utilized by GHRM. The errors labeled with the asterisks (*) are values corrected 
|47j from those in Ref. 





^^33(1232) 


Pii(1440) 


i?i3(1520) 


5ii(1535) 


5ii(1665) 


i^ir,(1680) 


i^37(1950) 


p (AJM) 
p (GHRM) 


-1.0±0.1 
-1.0±0.1 


-0.67 ±0.17 
(-0.62 ±0.16) 

'^•"^-0.20 


-0.84 ±0.17 

(-0.77 ±0.08) 

n 77+0.122 {,) 
U-' '-0.125 


-0.51 ±0.35 
-0.516«:?? 


-0.28 ±0.41 
-0.28t°:« (*) 


-0.27 ±0.08 
n 97+0.10 

^•^'-0.12 


-1±1 
-1±1 


n (AJM) 


-1.0±0.1 


-1.50 ±0.39 


-0.85 ±0.15 


-1.96 ±1.32 


-3.53 ±5.06 


-2.50 ±1.01 


-1±1 


d (AJM) 


-1.0±0.1 


-0.92 ±0.27 


-0.85 ±0.14 


-0.81 ±0.64 


-0.52 ±0.78 


-0.49 ±0.14 


-1±1 



using the ratio in Eq. (24), in analogy with Model II of GHRM |;14j. However, instead of 
fixing the parameters k^j so that the 77 and 7Z continuum pieces are equal [H] , we allow 
these to be determined by demanding that the 7Z structure functions be continuous across 
the boundaries of this region, that is, ai W = 3 GeV and Q^ = 2.5 GeV^. As we will see in 

rp J 

the following section, this places strong constraints on k^' , leading to significantly reduced 
uncertainties on the resulting value of 5Re D^^ . 

Finally, the 7Z structure functions in the DIS region, at W"^ > 4 GeV^ and Q^ > 
2.5 GeV^, are computed from the ABMll PDF parametrization ^011^ . The transformation 
from 77 to 'jZ is trivial at the parton level, amounting to a replacement of the quark electric 
charges e^ multiplying the universal PDFs by the weak vector charges gy. In the absence of 
'jZ structure function data at low Q^, the relative magnitude of the higher twist corrections 
to F2 was taken [48j to be the same as for F2"'. To account for this uncertainty, we therefore 
assign a conservative 5% uncertainty on F^ and F^ over the entire range of kinematics 
in Region HI. Since it is given by a difference of the F2 and F^ structure functions (see 



Eq. (14)), the longitudinal structure function F2 will necessarily have a larger relative 
uncertainty. 
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IV. PHENOMENOLOGICAL CONSTRAINTS 

As mentioned in the previous section, the central value of SfJeD^^ in Ref. [H] is given 
by the average of Models I and II, with the dominant nonresonant background contribution 
taken from Model II. If it were possible to reduce the background uncertainty, the error on 
the final 3fJe D^^ correction could also be lowered significantly. 

In their calculation of D^^ , GHRM [14j estimate the 7Z nonresonant background cross 
section by transforming the 77 cross section in the VMD+Regge model [29| according to 

"IT'' = I # 1 41". (29) 



0": 



77 



with the electromagnetic cross sections cfipf^ parametrized as in Eqs. (18a) and (18b), and 
the rescaling factor (c^tl/^'t'l) given by Eq. (24). The uncertainties on the 7Z cross section 



are obtained by comparing each Rj ratio in Eq. (24) with HERA data on exclusive vector 



meson electroproduction [IH] (c/. Fig. 13 of Ref. [E]), with the uncertainty taken to be the 
difference between the two. 



rp T 

The final contribution to the background error comes from the values of ^(j in Eq. (24) 



"c 

In the GHRM analysis |14| this term is equated with the electromagnetic continuum piece, 

assuming a 100% uncertainty. The resulting F2 structure function is illustrated in Fig. [6J 
as a function of both W"^ and Q^, and compared with the ABMll global fit [40]. Note that 
the uncertainty band on the GHRM VMD+Regge calculation includes only the continuum 
part of the background, and will be larger once the resonant uncertainty is included. The 
comparison clearly shows that the GHRM uncertainties are significantly larger than those 
typically obtained from global QCD analyses, especially in the region of intermediate W and 
Q'^ where both descriptions should be valid. Furthermore, as suggested already in Figs. [4] 
and |5j the central values lie systematically above the PDF parametrizations. 

Although the VMD model itself does not provide any additional constraints on the in- 
terference continuum contribution, we shall examine in this section the possibility of con- 

rp T 

straining K(j using existing knowledge of parton distributions, as well as recent data on 
parity- violating inelastic scattering from the Jefferson Lab E08-011 experiment [15]. These 
constraints will make it possible to reduce the overall uncertainty in 3fJe D^^ compared with 



those obtained in earlier analyses, Eq. (28). 
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FIG. 6: (color online) Comparison of the proton F2 structure function in the VMD+Regge 
model (Model II) of GHRM [14J (red dashed) with the ABMll global parametrization [lOj (green 
dotted), for fixed Q^ (top panels) and fixed W'^ (bottom panels). Note that the VMD+Regge 
model only includes uncertainties from the continuum part of the background, while the ABMll 
parametrization includes an overall 5% error. 

A. Constraints from PDFs 

In the deep-inelastic region at high W (W > 2 GeV) and Q"^ {Q"^ > 1 GeV^), structure 
functions can be described in terms of leading twist PDFs, with corrections from target mass 
and higher twist contributions included to account for residual, 1/Q^-suppressed nonpertur- 
bative effects. While a PDF-based description will eventually break down at low W and 
Q^, the region where the continuum contributions to the cross sections are relevant overlaps 
with the typical reach of global PDF parametrizations [JOl - l^ . One can therefore constrain 
the nonresonant part of the 7Z structure functions by requiring consistency of the model in 
the overlap region with the PDF parametrizations. 

Our fit of the parameters k^J involves equating the cross section ratios cf'^l/u'^l in 
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FIG. 7: (color online) Continuum parameters k^ (left) and k^ (right) fitted to the DIS data, 
parametrized by the ABMll global QCD fit [40], as a function of W"^ for fixed Q^ = 2.5 GeV^ 



r,L\ 



(red triangles), 6 GeV (blue squares), and 10 GeV (green circles). The average values {k^j ) are 
indicated by the solid lines, with the shaded band giving their uncertainty. Note that some of the 
points have been slightly offset for clarity. 



Eq. (24) with the structure function ratios computed from global QCD fits in the DIS 



region (see Eqs. (13) and (14)), 



jZ 



^77 



F, 



jZ 



F[ 



'77 



.7^ 



^77 



?lZ 



DIS 



;i77,7Z 



FT 



(30) 



DIS 



where the DIS structure functions F'^Y^ are taken from the ABMll parametrization |10] . 



As discussed in Sec. 



Ill 



T,L 



in fitting K(j in the DIS region, to be conservative we assume an 
overall 5% uncertainty on F^ , and a 40% uncertainty on Fl , which exceeds the uncer- 
tainties quoted in Ref. 0^ over the kinematics relevant for the 3fJe D^^ calculation. 

For the constrained fit we determine the values of k^! that minimize the x for each 
point in W^ and Q^, over a range of W^ values at fixed Q^ near the boundary between the 
DIS region (Region III) and the other regions in Fig. [2] To test the stability of the fitted 
Kj values with respect to the matching scale, we consider several different values of Q^ 
(Q"^ = 2.5, 6 and 10 GeV^). The resulting fits in Fig. ^indicate relatively mild dependence 
on the scale, which becomes negligible with increasing Q^ for k^, but with the expected 
larger uncertainties for Kq. 

The central values of Kq are computed by averaging over the three sets of Q^ values, 
and the uncertainty determined by taking into account both the W"^ dependence of the fits 



T,L 



and the PDF error. Because the kJ values at the different Q are correlated, performing a 
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simple x^ fit to all the sets may underestimate the errors. As a more reliable error estimate, 
we combine in quadrature the uncertainties arising from (i) the W"^ dependence, for which 
we take the average of the difference between the central values of the lowest and highest 
points for the Q^ set giving the strongest W'^ dependence (namely, for Q^ = 2.5 GeV^ for 
K^, and Q^ = 10 GeV^ for Kq); and (ii) the PDF error, the uncertainty for which is given 
by the data point with the largest error in the entire set (which occurs for Q^ = 2.5 GeV^ 
for both K^ and k^). The final fitted values of the continuum parameters are found to be 

K^ = 0.65 ±0.14, K^ = -1.3±1.7. (31) 



Compared with the uncertainties assumed by GHRM [15] our uncertainty on the transverse 
parameter k^ is about 5 times smaller, while that on the longitudinal parameter k,^ is almost 
two and a half times larger. However, the error on n^ has minimal effect on the 'jZ cross 
section at these kinematics because of the relatively small contribution of the longitudinal 
structure function. 

The resulting Fi] structure function with the constrained k^' values is shown in Fig. M 
for fixed Q^, ranging from Q^ = 0.05 GeV^ to 10 GeV^. The models of the 7Z structure 
functions are seen to match very well at the boundaries between the Regions I, II and III. 
As for the interference F2 structure function in Fig. pi only the continuum uncertainty 
is included in these examples; this allows a direct comparison with the uncertainty in the 
GHRM model input which dominates all other uncertainties. The comparison between 
Figs, y and p^ at the corresponding kinematics illustrates the significant reduction in the F2 
uncertainty that results from constraining the structure functions by the global QCD fits 
of PDFs. A similarly large reduction in the uncertainty can be seen in Fig. M for F2 as a 
function of Q^ at fixed W"^ values. 

The remaining uncertainty on the background contribution is associated with the -R^'^ 



-,T,L 



and R,' terms in Eq. (24). Following GHRM [H], we take the difference between these 
ratios calculated in the VMD+Regge model at Q^ = 7 GeV^ and the experimental vector 
meson cross sections from HERA [49j, assuming R'^ = R^ and R^ = R^ (see Fig. 13 of 
[13]). This uncertainty is then added in quadrature with the continuum uncertainty, along 
with the resonance contribution discussed in Sec. |III[ to obtain the total error on the 'yZ 
structure functions used in estimating 3fJe D^^ . 



The impact of the total uncertainty reduction is illustrated in Figs. [TO] and [TT] for the 
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FIG. 8: (color online) Proton F2 structure function versus W"^ at various fixed Q^ values for the 

low-W" CB fit [2i] (blue solid) and the high-W" VMD+Regge [29j (red dashed) and ABMll |iO] 

(green dotted) parametrizations. The boundaries between the Regions I, II and III are indicated 

by the vertical lines at VF^ = 4 and 9 GeV^. 
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FIG. 9: (color online) Proton F2 structure function versus Q^ at fixed W'^ = 4, 6, 9 and 12 GeV^ 
for the CB fit [24] (blue sofid), the ABM 11 PDF parametrization [40j (green dotted), and the 
VMD+Regge model [29] (red dashed), with the boundaries between Regions I, II and III indicated 
by the vertical lines at fixed Q^. 



parity-violating inelastic asymmetry for the proton, 
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where y = v/E is the fractional energy transferred to the target. In addition to the vector 
F12 structure functions, the asymmetry A-p\ depends also on the axial-vector F^ structure 
function. For the the resonance contribution to F^ we use the parametrization of the axial- 
vector transition form factors of Lalakulich et al. [5IH5S]- For the background we follow 
Ref. [31] and rescale the electromagnetic cross sections [2^ by the average of the x — )■ 
and SU(6) quark model limits, which gives F^ = 5/3 F^^. (Note that for the deuteron this 
average becomes F^ =9/5F^^.) 



The asymmetries calculated in the AJM and GHRM models are shown in Fig. 10 at an 
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FIG. 10: (color online) Proton parity-violating inelastic asymmetry ^pv/Q^ as a function of W ^ 
at fixed incident energy E = 0.69 GeV and Q^ = o.34 GeV^, for the GHRM Model II [14J (left) 
and the AJM model (right). The data point at W = 1.18 GeV (black circle) is from the Jefferson 
Lab GO experiment |50| . 
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FIG. 11: (color online) Proton parity-violating inelastic asymmetry Apy/Q^ as a function of W , 
at fixed incident energy E = Q GeV and Q^ = 2.5 GeV^, for the GHRM Model II [H] (left) and 
the AJM model (right). The asymmetry computed directly from PDFs |40j is represented by the 
green band. 

incident energy E = 0.69 GeV and Q^ = 0.34 GeV^, corresponding to the kinematics of 
the recent GO measurement at Jefferson Lab near the A resonance region [30]. The central 
values of both models agree well with the data, although the experimental uncertainty is 
too large to enable meaningful constraints to be placed on the jZ structure functions. The 
constraint on the Kq value from matching to the DIS structure functions in the AJM model 
renders the uncertainty band somewhat smaller than the GHRM uncertainty [H] at higher 
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values of W. (Note that the uncertainty on Apy is computed by taking the upper and lower 
values of the input 'jZ structure functions, and is therefore asymmetric.) 

The difference in the error bands becomes more pronounced at larger Q^, as seen in 



Fig. 11 at i^^ = 6 GeV and Q^ = 2.5 GeV^, which are representative of typical kinematics at 



Jefferson Lab (see Sec. IV B below). Here the uncertainty on the GHRM model asymmetry at 



W ^ 2 GeV is around 4 times larger than the corresponding uncertainty on the constrained 



A JM model asymmetry. For comparison, we also show in Fig. [TT] the asymmetry computed 
directly from PDFs [ID] in the region W > 2 GeV where a partonic description is expected 
to be valid. 

The uncertainty in the PDF-based calculation is slightly smaller than, but qualitatively 
similar to, that in the AJM model, while the GHRM model uncertainty is significantly 
overestimated in the region of overlap. We stress that although the DIS region makes only a 
modest contribution to 3f?e D^^ , the requirement that the ■yZ cross sections match across the 
DIS-resonance region boundary imposes strong constraints on the 'jZ structure functions 
also at lower W and Q^. In the following section we confront this against new data on 
parity-violating electron-deuteron scattering in the resonance region. 



B. Deuteron asymmetry 



The E08-011 experiment |T5] at Jefferson Lab recently measured the parity- violating 
asymmetry in inclusive electron-deuteron scattering over a range of W and Q^ in both the 
resonance and DIS regions. While the DIS region data are currently still being analyzed 
[5^ . the available resonance region data [55] can be used to provide an independent test of 
the procedure for estimating the 'jZ structure functions. This is particularly important for 



3fJe Di;^ , since the integrals in Eq. ( 12 ) are dominated by Region I in Fig. 2 

The measured parity-violating asymmetry Apy, scaled by l/Q"^, is shown in Fig. 12 at 
W = 1.26, 1.59, 1.86 and 1.98 GeV, with Q^ values ranging from 0.76 GeV^ to 1.47 GeV^. 
(The l/Q"^ scaling factor enables the various points to be shown on the same graph.) The 
deuteron asymmetries in the AJM model are computed with the continuum parameters 
constrained by the DIS region structure functions computed from global PDFs [40J, as for 



the proton asymmetry in the previous section (see Fig. 11). The resulting fit gives for the 
transverse continuum parameter Kc<{d) = 0.79 ±0.05, and is in excellent agreement with the 
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FIG. 12: (color online) Deuteron parity-violating asymmetry Apy/Q"^ as a function of W for 
incident electron energy E = 4.9 GeV (left) and E = 6.1 GeV (right). The data points from the 
Jefferson Lab E08-011 experiment flS] at W = 1.26 (green square), 1.59 (red circle), 1.86 (blue 
triangle) and 1.98 GeV (black diamond) correspond to average values of Q^ = 0.95, 0.83, 0.76 and 
1.47 GeV^, respectively. The AJM model uncertainties (inner dashed band) are constrained by 
matching the continuum parameters n^j (d) to the DIS region jZ structure functions [40j, and 
are compared with those computed with errors on k^j (d) of 100% (outer dotted bands) and 25% 
(inner dotted bands). 



E08-011 data [15] for all kinematics, except at the A region point at Q^ = 0.95 GeV^, where 



it lies slightly below the data. Since the calculation of the A resonance contribution to A 



PV 



relies only on isospin symmetry and the conservation of the vector current, its uncertainty is 
smaller than that for higher-mass resonances. The discrepancy may reflect stronger isospin 
dependence of the nonresonant background for A production |56], although the difference 



is at the < 2a level. Also, as seen in Fig. 10 above, the models agree well with the GO data 



[50] in the A region, albeit within larger errors. 

By using the longitudinal structure function from the global QCD fit in Ref. [10], we find 
for the longitudinal continuum parameter Hc{d) = 0.2 ± 3.4. Although the specific imple- 
mentation of the CB parametrization [24j prevents this uncertainty from being propagated 
directly into Apy, we nevertheless can use the Kfj values for the proton to ensure that 
the uncertainty in the longitudinal piece is taken into account. For comparison, we also 



show in Fig. 12 the uncertainty that would be obtained with a similar 100% error on the 
continuum parameters as was assumed by GHRM for the proton, with the VMD+Regge 
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FIG. 13: (color online) As in Fig. 12, but with the AJM model asymmetries (solid) and their 
uncertainties (dashed) constrained by the E08-011 data [55j. Note the different scale on the 
ordinate to that in Fig. [12] 



model [2n] used for the entire kinematic region |14j. In this case the uncertainties on Apy 
in the ly > 1.8 GeV region are ~ 6 times larger than the AJM model asymmetries. Using a 
reduced 25% uncertainty on K,c{d) results in asymmetries with a significantly smaller error 
band, which is nevertheless slightly larger than in the AJM model. 

As a check, the parameter K^{d) was also constrained by performing a x^ fit to the E08- 
011 data points. This fit constrains the dominant, transverse continuum parameter to be 
K^{d) = 0.69 ± 0.13. (Omitting the A datum from the fit would yield a marginally larger 
value, Kc{d) = 0.72 ± 0.13.) For the longitudinal contribution, the CB parametrization 
of the deuteron structure function provides only -F^'^, while F2'^ is obtained through the 



longitudinal to transverse cross section ratio a]^ /a'^ (see Eq. (21)), with the deuteron ratio 



assumed to be the same as for the proton. Within this parametrization, a direct constraint on 
Kcid) £^s for the proton case is therefore not possible. However, as for the PDF-constrained 
asymmetry, we can still propagate the uncertainty on ai/aT through the final asymmetry 



T,L 



by including the uncertainties in the K(j values of the proton which serve as inputs into the 

af/af 



ratio. 
The resulting asymmetries are again in very good agreement with the E08-011 data, as 



is seen in Fig. 13 Moreover, the uncertainties (dashed curves) are 3 — 4 times smaller 
in the W^ > 1.8 GeV region than those obtained by assuming a 100% uncertainty on the 
parameters, and remain smaller than even for the reduced, 25% uncertainty case. The 
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TABLE II: Parity- violating deuteron asymmetries in the AJM model at the kinematics of the 
E08-011 experiment [15j. The asymmetries are computed with the continuum parameters Kq (d) 
constrained by the E08-011 data, or by matching to the DIS region described in terms of PDFs. 
Note that the points marked with asterisks (*) are predictions. 



E W Q^ 

(GeV) (GeV) (GeV^) 



Apy/Q^ (ppm GeV- 



PDF constraint 



E08-011 constraint 



4.9 
4.9 
4.9 
6.1 



1.26 
1.59 
1.86 
1.98 



0.95 
0.83 
0.76 

1.47 




-93.1 



+8.8 
-9.0 



-80.1 

-82.4 



-flO.l 



10.3 

r7.9 

-8.0 



-79.2+8| 



6.1 
6.1 
6.1 



2.03 
2.07 
2.33 



1.28 
1.09 
1.90 



-84.9+^:^ (*) 
-85.2+6;2 (*) 
_g2.7+6.3 (*) 



-70 7+8-4 (*) 
'^•'-8.6 

-80.3+1^ (*) 
-76.5:^9.3 (*) 



consistency between the data and the results given by the constrained expressions gives us 
confidence in the reliability of the 7Z structure functions in the AJM model in the region 
of low to intermediate W and Q^ that is of greatest importance for the 3fJe Dl^^ calculation. 
Finally, the values of the calculated asymmetries and their uncertainties, using both the 
resonance region data and the PDF constraints, are summarized in Table [TT] at each of 
the kinematic points from the E08-011 experiment [15]. In addition, we list the AJM model 
predictions for Apy at the measured DIS region points a.tW>2 GeV (marked by asterisks), 
which will be discussed further in the next section. 



V. RESULTS 



A. jZ box corrections for Qweak 



The detailed examination of the 'jZ interference structure functions and their uncertain- 
ties, constrained by data in the DIS region and parity-violating asymmetries in the resonance 

J7Z 



region, allows us to compute the '^mnY^ correction in Eq. (12), and through the dispersion 
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FIG. 14: (color online) Energy dependence of the contributions to KeDY^ from the various regions 
in W and Q^ displayed in Fig. [2] in the AJM model (top), and the breakdown of Region I into its 
resonant and nonresonant background components (bottom). 



relation (pi) to make a reliable determination of the 7Z box correction to Q^- The depen- 
dence of 3f?e n^^ on the incident energy E is illustrated in Fig. 14l which also shows the 
individual contributions of the various W and Q^ regions in Fig. ^ 



At low energy {E <\ GeV), the total correction SfJeD^^ is dominated by the low-iy, 
low-Q^ region (Region I in Fig. ^. As found in earlier analyses (UHlll [22], the resonant 
contribution (mainly from the A(1232) resonance) peaks at around E ^ 0.7 GeV, and grad- 
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TABLE III: Contributions to 3f?e D^^ from various regions in W and Q in the AJM model (see 
Fig. [21) at the Q^eak energy E = 1.165 GeV. 



Region 


Ken:^^ (xlO-3) 


I (res) 


2.18 ±0.29 


I (bgd) 


2.46 ±0.20 


I (total) 


4.64 ±0.35 


II 


0.62 ±0.05 


III 


0.35 ±0.02 


total 


5.60 ± 0.36 



ually decreases at higher energies. The nonresonant and resonant components of Region I 
are approximately equal at i? ~ 1 GeV, with the nonresonant part growing with increasing 
energy. The higher-VT , higher-Q^ regions play a relatively minor role in the dY^ correction, 
with Regions II and III contributing ^ 20% and 10% of the total, at E = 3 GeV, respectively. 
At the Qweak energy, E = 1.165 GeV, the breakdown of the ^ReD^^ correction into its 



individual contributions is summarized in Table pTT| Including uncertainties from all regions, 
the total correction is found to be 



3fJen;^ = (5.60 ±0.22 [bgd] ±0.29 [,es] ±0.02 [Dis]) x 10" 



(33) 



where the uncertainties listed are from the nonresonant background, the resonances, and the 
DIS region, respectively. Adding the errors in quadrature gives 3fJe D^^ = (5.60±0.36) x 10~^ 
at the Qweak energy. The ^ 7% relative uncertainty on this correction remains largely energy 
independent, even at large energies, where the contributions from larger W and Q^ become 
more important; since the structure functions are constrained by DIS data, the uncertainty 
in 3fJe D^^ does not grow with E. 

The AJM model value of the jZ box correction is similar to the result, 3fJe D^^ = (5.40 ± 
0.54) X 10~^, obtained using the 7Z structure functions from Region II extended over all 
kinematics, as in the GHRM Model II [H], but with the kJ parameters constrained by 
matching to the DIS region structure functions [40j . This constraint renders the uncertainty 
~ 4 times smaller than that in Ref. [H], but still slightly larger than in the AJM model 
calculation. 
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FIG. 15: (color online) Predictions for the parity-violating deuteron asymmetry Ap^/Q'^ as a func- 
tion of W (solid) for the DIS region kinematics of the Jefferson Lab E08-011 experiment [15] at 
Q2 = 1.28 GeV2 (green), 1.09 GeV^ (red) and 1.90 GeV^ (blue) (see also Table|lll). The uncertain- 
ties (dashed) are computed in the AJM model with the continuum parameters Kq constrained by 
DIS structure functions (left), and by the E08-011 resonance region data (right). The predictions 
at the experimental W values [15] are shown as pseudo-data points (open symbols). 



B. Predictions for parity- violating asymmetries 

The 7Z structure functions can be further constrained by additional parity-violating 
asymmetry data from the E08-011 experiment at Jefferson Lab [I5]. The deep-inelastic 
region data are currently being analyzed [51], and the predictions from the AJM model are 



shown in Fig. 15 as a function of W for the three experimental Q^ values (see also Table ^ 



The uncertainties on the predictions are computed both by fitting the continuum parameters 

T' T • i ■ 'i . n . ^^ 

." ' -f-/-* r \r\ r-\ lll^ c^r -vi -i ftr i -i -v ^-\ tit -irt f-ir i /~\-ir\ i~t I /I I 1 1 n -tn /~\ -t-|^/~\ UlllW llTn v/~\i~ij— vv-i(-i -irt f-i/-\ ■v/-\rfi /"w^t /~\ i-\r i~\ ITl^l I MO 



K(j to the DIS structure functions |40j and the E08-011 resonance region data 
asymmetries with the E08-011 data constraints are marginally higher than those with the 
parameters constrained by PDFs, with slightly larger uncertainties. As for the resonance 



region comparison in Figs. [12] and [13} these uncertainties are f« 4 — 5 times smaller than they 
would be without the constraints on K^j , assuming 100% errors along the lines of the proton 
calculation in Ref. [H]. The upcoming data will therefore be extremely useful in determining 



the uncertainties on the 7Z structure functions and on the resulting 3fJe D^^ correction. 



A further constraint will be provided by the inelastic Qweak measurement [IB], which 
was a special run of the Qweak experiment tuned to the inelastic region at an average W = 
2.23 GeV. The AJM model prediction for the proton asymmetry Apy and its uncertainty are 
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FIG. 16: (color online) AJM model prediction for the proton parity-violating asymmetry ^py as 
a function of W for the Qweak inelastic measurement [16j at Q^ = 0.09 GeV^ (solid line and open 
symbol). The AJM model uncertainties (dashed) are compared with those from the GHRM model 
with 100% uncertainty on the continuum parameters (dotted). 



shown in Fig. 16 , where we find Ap 



PV 



-7.8±0.6) ppm at the experimental Q^ = 0.09 GeV^ 



T,L 



value. The uncertainty in the AJM model, with the continuum parameters n^j constrained 
by the DIS structure functions, is ~ 2 times smaller at the inelastic Qweak kinematic point 
than that from the GHRM model [H] without these constraints. Note also that in the 
resonance region, W ^ 1.5 GeV, the uncertainty in the GHRM model almost doubles 
by taking extrema values instead of the more conventional addition in quadrature. The 
inelastic Qweak, and similar measurements of the parity-violating inelastic asymmetries, will 
be valuable for constraining the 7Z structure functions and the 3fJe D^^ corrections in the 
future. 



VI. CONCLUSION 

We have performed a comprehensive analysis of the 7Z box contribution to the forward 
electron-proton elastic parity-violating asymmetry. Our primary result is a new determi- 
nation of the uncertainty on ^ed^^ at the beam energy of the Qweak experiment. In 
comparison with previous estimates, we report a significant reduction in this uncertainty, 
driven largely by data on structure functions in the DIS region, and measurements of parity- 
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violating asymmetries in the resonance region. 

To isolate the dependence on the various inputs required in the evaluation of ^e D^^ 5 '^^ 
have divided the dispersion integral into three kinematic regions. Region I, which includes 
resonance contributions at low W and Q^, is identified to totally dominate the value of 
3fJe n^z ■ The total uncertainty is therefore largely driven by how well the 7Z interference 
structure functions F^ can be constrained in this region. 

The resonance region 7Z structure functions are determined by an isospin transformation 
of the corresponding 77 structure functions. The input F^'^ functions are determined by 
a fit [24J to the world's inclusive electron-nucleon scattering data in terms of resonance 
contributions and a nonresonant background. For the resonance components, the isospin 
transformation can be performed using the conservation of the vector current and the isospin 
dependence of the couplings, as reported by the PDG, with relatively modest contribution 
to the overall uncertainty. For the background, following the approach of Ref. [H], the 
transformation is estimated using a prescription based on the VMD model [22] • For the 
low-mass vector meson components the isospin rotation is determined by isospin symmetry 
of the electroweak interactions, while the transformation of the high-mass continuum part 
is not fixed within the VMD formalism, and consequently contributes a larger uncertainty. 

At larger Q^ values {Q"^ > 1.5 GeV^) the continuum piece totally dominates the nonres- 
onant background. We use this fact to constrain the continuum component of the isospin 
rotation by matching this to the DIS structure functions in the transition region. The model 
dependence from using a particular continuum form at lower Q^ (away from the PDF con- 
straint) is less important, since this region is dominated by the low-mass vector mesons 
p, u and 0. It is the constraint on this rotation that drives the significant reduction in 
uncertainty in the present AJM model as compared to that reported by GHRM [Ti] . 

Combined with the relatively well-determined contributions from Regions II and III at 
higher W and Q^ (see Fig. p|, we find the final value for 7Z correction to be SfJeD}^^ = 
(5.60 ± 0.36) X 10"'^. Importantly, this precision maintains confidence in the interpretation 
of the Qweak experiment as a standard model test. 

The reliability of our constraint procedure has been confirmed by a comparison with 
the corresponding inclusive 7Z interference asymmetries recently measured on the deuteron 
by the E08-011 experiment at Jefferson Lab [23]. Conversely, using the E08-011 resonance 
region data as a constraint on the 7Z structure functions, the resulting asymmetries are 
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found to be very similar to those in the AJM model with the PDF constraints, albeit with 
slightly larger uncertainties. Upcoming data on the deuteron asymmetry in the DIS region 
^^ should reduce these uncertainties. 

Beyond this, the most promising means by which one could further constrain the 'yZ 
structure functions would be to perform a systematic experimental study of parity- violating 
electron scattering on hydrogen across Region I. While the recent deuterium measurements 
[55] have proven useful in providing confidence in the procedure of matching to PDFs at 
intermediate Q^ and W, because the deuteron requires a knowledge of the neutron structure 
function as well as of the proton, this has limited value as a means to reduce the uncertainty 
in F^ . A dedicated study of the proton itself would directly constrain the model and lead 
to a reduction in the uncertainty of the radiative correction arising from the 7Z box. 
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